Transport of interacting electrons in arrays of quantum dots and diffusive wires 



O 

o 

(N 



Dmitri S. Golubev^''^ and Andrei D. Zaikin^'^ 
^Institut fur Theoretische Festkorperphysik, Universitat Karlsruhe, 76128 Karlsruhe, Germany 
Forschungszentrum Karlsruhe, Institut fur Nanotechnologie, 76021, Karlsruhe, Germany 
^I.E. Tamm Department of Theoretical Physics, P.N. Lebedev Physics Institute, 119991 Moscow, Russia 

We develop a detailed theoretical investigation of the effect of Coulomb interaction on electron 
transport in arrays of chaotic quantum dots and diffusive metallic wires. Employing the real time 
path integral technique we formulate a new Langevin-type of approach which exploits a direct 
relation between shot noise and interaction effects in mesoscopic conductors. With the aid of this 
approach we establish a general expression for the Fano factor of ID quantum dot arrays and derive 
a complete formula for the interaction correction to the current which embraces all perturbative 
results previously obtained for various quasi-OD and quasi- ID disordered conductors and extends 
these results to yet unexplored regimes. 
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I. INTRODUCTION 

Recently a profound relation was elucidatediiSiSi^ 
between full counting statistics (FCS)^ and electron- 
electron interaction effects in coherent mesoscopic con- 
ductors. In particular, it was demonstrated^ that the 
leading interaction correction to the current through such 
conductors is determined by the second cumulant of the 
current operator 52, i.e. by the power spectrum of the 
shot noise^. The interaction correction to the shot noise 
52 was in turn found^ proportional to the third cumulant 
of the current operator ^3 . Even more generally, one can 
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sho-i 



that the lowest order interaction correction to the 



n-th current cumulant 5„ is determined by Sn+i for all 
values of n. Since the current cumulants in the absence 
of interactions can be conveniently analyzed within the 
FCS formalism^, the above observations provide a great 
deal of information about the effect of electron-electron 
interactions as well. 

In order to investigate the influence of interactions on 
higher current cumulants it is in general necessary to em- 
ploy a complete expression for the effective action of a co- 
herent scatteresSi^. At the same time the resultsi for the 
first cumulant, i.e. the relation between the leading in- 
teraction correction to the current and the shot noise can 
be understood already within a simple and transparent 
theoretical framework of quasiclassical Langevin equa- 
tions. In the case of a single coherent scatterer shunted 
by some linear Ohmic resistor Rs these equations take a 
remarkably simple form 
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Here C is the scatterer capacitance, (p/e — V is the fluc- 
tuating voltage across the scatterer and Vx is the total 
voltage applied to the system "scatterer -|-shunt" . As usu- 
ally, one describes the scatterer by a set of conducting 
channels with transmissions Tk. The scatterer conduc- 
tance is then defined by means of the standard Landauer 
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^(t) is the noise of the scatterer, characterized by the 
correlator 
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is the Fano factor, and fs(0 the equilibrium noise of 
the shunt with the correlator 
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The whole approach based on Eqs. (|1I5|I is applicable 
either at sufficiently high energies or, more importantly, 
if at least one of the two dimensionless conductances, 
g = 27r/e^i? and/or gs = 2-K/e^Rs, remains large. The 
above Langevin equations make the relation between the 
interaction correction to the current and the shot noise 
completely transparent demonstrating that the former 
can be derived just if one accounts for the noise contri- 
bution in the equation describing the balance of currents 
across the scatterer. 

The above simple approach may hold only for rela- 
tively compact scatterers, in which case the description 
of interaction effects with the aid of the uniform in space 
fluctuating field is sufficient. For spatially extended 
conductors the coordinate dependence of this field cannot 
anymore be disregarded and the whole analysis needs to 
be modified. This modification is trivial if one considers 
an array of scatterers connected between each other by 
relatively big metallic grains. Assuming that the electron 
distribution function in each of these grains remains in 
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equilibrium one can describe the n-th scatterer by its own 
fluctuating field ipn which obeys the same set of Langevin 
equations Ql, ©■ For arrays of tunnel junctions this ap- 
proach was employed in Ref. 0. The corresponding gen- 
eralization of the results^ to the case of arbitrary scat- 
terers just requires modification of the Fano factor in the 
noise correlator J^Jl. 

The condition of local equilibrium inside metallic 
grains implies that the corresponding electron dwell time 
td between two adjacent scatterers should be much 
longer than the inelastic relaxation time Tin. If this con- 
dition is violated, the electron distribution function is 
driven out of equilibrium and the whole consideration be- 
comes more complicated. In the case of a quantum dot 
formed by two arbitrary scatterers the latter situation 
was analyzed in details in Ref. |^ and also in Refs. ^|9j for 
the case of chaotic dots. In all these works it was demon- 
strated that in the limit of large conductances and in the 
voltage biased regime the interaction correction to the 
conductance tends to saturate at temperatures/ voltages 
below 1/td- This implies that for finite values of td 
highly conducting quantum dots should show metallic 
behavior down to zero temperature. 

It is important to emphasize that this observation 
holds only provided the voltage source is attached di- 
rectly to the quantum dot, i.e. the external impedance 
is equal to zero. For non-zero external impedances volt- 
age fluctuations lift the conductance saturation, and the 
amplitude of the interaction correction keeps increas- 
ing with decreasing T even at temperatures well be- 
low 1/t£). In this regime the interaction correction was 
found* to scale linearly with the total Fano factor of the 
quantum dot and to depend logarithmically on tempera- 
ture/voltage for sufficiently large external impedances or 
if this impedance is purely Ohmic. 

For similar reasons no saturation of the interaction cor- 
rection at energies below 1/td should be expected for 
chains and arrays of quantum dots. Recently this situa- 
tion was analyzed diagrammatically^" in the case of gran- 
ular tunnel junction arrays. Indeed, it was found that the 
interaction correction increases with decreasing temper- 
ature both above and below the inverse dwell time in 
individual grains. MT < 1/td the authorsiS recovered 
exactly the same expression for the interaction correc- 
tion as that known in the case of diffusive conductorsii. 
This equivalence is by no means surprising if one bears 
in mind the fundamental relation between the interac- 
tion correction and the shot noise on one hand, and the 
resultsiSii^ on the other hand, which demonstrate that 
the shot noise of a sufficiently long array of tunnel junc- 
tions is equivalent to that of a diffusive wire. Extending 
these arguments to arbitrary scatterers, with the aid of 
the resultsJS^ one can anticipate that at sufficiently low 
energies (< I/tu) and large spatial scales the interac- 
tion correction should be described universally for any 
array of quantum dots and ultimately for any mesoscopic 
conductor in the metallic regime. This universality will 
indeed be demonstrated below. 



The main goal of the present paper is to generalize 
the simple Langevin equation approacl^ to situations 
in which relaxation of the electron distribution func- 
tion occurs at much longer time scales as compared to 
the electron dwell time between two adjacent scatterers 
Td ^ Tin. Although the distribution function may sig- 
nificantly deviate from the Fermi function, it is possible 
to account for these deviations within the (generalized) 
Langevin equation analysis and to formulate a closed set 
of equations which fully determine the interaction cor- 
rection to the I — V curve of disordered conductors. 

The structure of the paper is as follows. In Sec. II we 
will specify the model of a disordered metallic conduc- 
tor and present a phenomenological derivation of the ba- 
sic Langevin equations for our problem. This derivation 
will be carried out with the aid of simple and transpar- 
ent physical arguments which make the whole approach 
easy to understand without going into technical details. 
A more advanced analysis employing the effective ac- 
tion technique will be described in Sec. III. This anal- 
ysis provides rigorous justification for our phenomeno- 
logical derivation and allows to illustrate a useful rela- 
tion between our technique and the classical Boltzmann- 
Langevin approach^. In Sec. IV we will probe our 
Langevin technique by explicitly deriving the shot noise 
spectrum and the Fano factor for arrays of chaotic quan- 
tum dots in the absence of interactions. The remain- 
der of the paper will be devoted to the analysis of the 
leading interaction correction to the current in arrays of 
quantum dots and mesoscopic diffusive wires. In Sec. 
V we will derive the general expression for this correc- 
tion which then will be applied to homogeneous arrays 
of quantum dots in Sec. VI. In the latter case we will es- 
tablish a complete analytic form of the interaction correc- 
tion and present the corresponding simplified expressions 
in a number of important limits. Our general formula, 
Eq. (|63() . embraces all previous resultgi^iSi^iifliiiJ'* ob- 
tained in various types of quasi-OD and quasi-lD disor- 
dered conductors, allows to establish a transparent re- 
lation between these results and to extend them to yet 
unexplored regimes. A brief analysis of an additional ef- 
fect of external leads will be presented in Sec. VII. We 
will then discuss our results and conclude the paper in 
Sec. VIII. 



II. THE MODEL AND PHENOMENOLOGICAL 
ANALYSIS 

We shall consider a chain of — 1 quantum dots as 
it is shown in Fig. 1. Each dot can be viewed as an 
island in-between two scatterers/barriers which in turn 
connect adjacent quantum dots. Electrons can enter the 
dot through one of the barriers, spend some time there 
propagating between the barriers, possibly being scat- 
tered at the barriers, outer walls or otherwise, and finally 
leave the dot through another barrier. In what follows 
we will adopt the model of chaotic quantum dots. 
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FIG. 1: ID array of chaotic quantum dots. The array consists 
of A'' — 1 dots and A'' barriers. The n-th dot is characterized 
by mean level spacing 5n and gate capacitance Cgn- The n- 
th barrier is described by its Landauer conductance 1/Rn, 
capacitance C„ and Fano factor f3n. The array is placed in- 
between two big metallic reservoirs which are connected to 
the voltage source via Ohmic resistor Rs- 



Each of the N barriers will be described by its Lan- 
dauer conductance 1/Rn ~ (g^/tt) X]fc ^i"^ capacitance 



C„ and Fano factor /3„ - E;c?^i"^(l " T^rO/Efc?^,! 

(n) 

where is the transmission of the k-th conducting 

mode in the n-th barrier. We also define dimcnsionless 
conductances of the scatterers gn = 27r/e^i?„. In what 
follows we will assume that each scatterer has many con- 
ducting channels and that its dimcnsionless conductance 
is large (?„ ^ 1. The n-th dot will be characterized by 
the mean level spacing Sn = I/NqVu, where Vn is the dot 
volume and A^o is the density of states at the Fermi level. 
For the sake of generality we will also assume that each 
dot has an additional capacitance to the ground Cgn ■ Fi- 
nally, the first and the last scatterers are connected to 
two big metallic reservoirs which in turn are connected 
to the voltage source via external leads with an Ohmic 
resistance Rs- 

An important assumption concerns the spatial depen- 
dence of fluctuating voltages in our system. Similarly to 
Ref. 's' we will allow for voltage drops Vn{t) only across 
the barriers, while inside the dots voltages can depend ar- 
bitrarily on time but not on the spatial coordinates. This 
assumption is usually well satisfied for metallic dots con- 
sidered here. In the leads the voltage fields are assumed 
to vary slowly in space. In the course of our derivation 
we will essentially neglect all mechanisms of inelastic re- 
laxation which are not related to electron-electron inter- 
actions. We will also disregard weak localization effects 
which can be easily suppressed, e.g., by externally ap- 
plied magnetic field. 

We will now proceed with our phcnomcnological anal- 
ysis of the above model. 



also assume that the electron energies are measured with 
respect to the unique reference energy for the whole array. 
In this case the noise of the n— th scatterer 5„(ii,t2) = 
{^7i{ti)(.n{t2)) takes the form^ 



1 f duj 



2i?„ J 2tt 

x{Pn[fn-l{E)hn{E + OJ) + fn{E + Lj)hr,-l{E) 

+fn{E)h„-i{E + Lj) + f„-i{E + uj)hn{E)] 

+ - Pn)[fn-liE)hn^l{E + LU) 
+fn-l{E + L0)K-l{E) + fn{E)K{E + LJ) 

+UE + Lu)KiE)]}, (6) 
where hn{E) = 1 — fn{E). Let us define the function 
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dE 
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-lE{ti-t2) 



(l-2/„(i?)). (7) 



In stationary situations this function depends only on 
the time difference ti — 12 and it also obeys the condition 
Gn{t2,ti) = C?*(fi,t2)- In equilibrium the distribution 
function equals to the Fermi function 



fF{E) 



1 



exp(f ) 
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Substituting this function into Eq. ((Jj) one finds 
Gn{ti,t2) = -iTRe [l/sinh7rr(ii - ^2 + «e)] , where e is 
an infinitesimal positive constant. At ti — )■ t2 one gets 
Gn{ti,t2) — iRe [l/7r(ti — ^2 + ie)]- This analytical 
property turns out to be general, i.e. it equally applies 
both to equilibrium and non-equilibrium situations. Ex- 
pressing Eq. via Gn{ti,t2), we obtain 



Sn{ti, ^2) 



7ri?„((ii-i2)2+e2)2 



(9) 



2Rn 



t2)G'n(^2i ^l) + Gn{ti,t2)Gn-l{t2,ti)^ 



^(l-/3n) 
2Rn 



[\Gn-l{tl,t2)\^ + \Gn{tl,t2)\' 



Although the formula ^ has been derived under the as- 
sumption that the distribution function does not depend 
on time, we will show later that it remains valid also in 
non-stationary situations. In the latter case the function 
Gn{ti,t2) can be understood as the Keldysh component 
of the quasiclassical Usadel Green function. 



A. Noise correlator 

As a first step we will specify the general expression 
for the noise correlator needed for our derivation. Let us 
assume that the electron distribution function fn{E) in 
the ri-th dot does not depend on time for all n but may 
deviate arbitrarily from the Fermi function. Below we 



B. Kinetic equation 

Our next step is to derive the kinetic equation for the 
function G'„(ii, 12). For simplicity we again start from the 
stationary situation, in which case fluctuations of volt- 
ages in our system can be neglected. In what follows we 
will assume that both /„(£') and G'„(ti, t2) do not depend 
on coordinates inside the n— th dot. The total number of 
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electrons with energies in the interval [E, E + dE] in the 
n-th dot is 2NQVnfn{E)dE , where the factor 2 accounts 
for spin degeneracy. This number may change in time 
only if electrons leave and/or enter the dot trough the 
left (n— th) and the right (n + 1— th) barriers. One finds 



2dEdUt,E) 

Sn dt 



TV. 



(10) 



The escape rate from the n-th dot and the transition 
rate to this dot - both through the n— th scatterer - are 
respectively 

r„-i,,i = c/„(5„/47r, Tn,n-i = g„(5„-i/47r. (11) 
Then for J°"* one obtains 

jr* = 2(r„+i,„ + r„_i,„)/„(i;)d£;/5„, 

and similarly for J™. Combining the above simple ex- 
pressions we arrive at the kinetic equation 



5„ dt 



-{g-n + 9n+l)fn + 9nfn-l + (12) 



where 1 < n < iV — 1. The boundary conditions to 
this kinetic equation are set by the requirement that the 
distribution functions in the left and the right reservoirs, 
/o and /tv, are equal to the Fermi function, i.e. 



fo (E) = fp (E) , fN (E) ^fpiE- eV) . 



(13) 



Here and below V is the total voltage applied to the array. 
We note that the function 1 — 2/„ also satisfies Eq. (I12II . 

The kinetic equation for the function G'„(ti, t2) can be 
obtained from Eq. ((T^ if we identify < = (ti + i2)/2, 
introduce s = ti — t2 and make the Fourier transform 
of lO by taking the integral / |f e-^^''(. . . ). Then we 
obtain 



= -i9n+ gn+l)Gn{t,s) + gnGn-l{t, s) 



+ gn+lGn+lit, s) 



(14) 



As we have already pointed out Eq. H14|l applies only 
in stationary situations. A proper generalization of this 
equation for non-stationary cases can be achieved with 
the aid of general gauge invariance arguments which yield 

= + 5'n+l)G„(il,<2) + 5nG,i-l(tl, ^2) 

+ 9n+lGn+l{tl,t2), (15) 

where we defined <i>„(t) = J2^=i^^ji^)- This kinetic 
equation holds for arbitrary time dependent voltages. As 
before, the boundary conditions to this equation read 



Go(ii,i2) = -iTRe 



1 



sinh7rr(ti — t2 + «e) ' 



GN{ti,t2) = -^^e-'"^(*l-*^)Re 



1 



sinh7rT(ii — t2 + ie) 
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C. Balancing fluctuating charges and voltages 

In order to complete our simple analysis we formu- 
late the standard circuit theory equations, which allow 
to include the effect of charge accumulation in quantum 
dots. Let us define the fluctuating excess charge in the 
n—th dot g„. We assume that all quantum dots are well 
described by the capacitance model, in which case one 
finds 



Gn+lVn+l — CnVn 



(17) 



Here Cn is the capacitance of the n—th barrier and Ggn 
is the capacitance of the n—th dot to the ground. The 
current /„ is formed by the sum of three different terms, 
namely the standard Ohmic term Vn/Rn, the "kinetic" 
term r„^„_ig„_i - r„_i^„g„ and the noise term 



Rn 



(18) 



The variation of the charge qn is in turn determined by 
the currents flowing through the n—th and n + 1— th bar- 
riers. We get 



Qn — ^71 ^n-\ 



(19) 



Finally, the sum of all fluctuating voltages Vn should be 
equal to the total applied voltage. 



N 



J2Vn = V. 



(20) 



Eqs. (|9I15I17I20|I form a complete set of equations, 
which allow to find the first order interaction correction 
to the I — V characteristics for an array of quantum 
dots. These equations represent a straightforward gen- 
eralization of the Langevin approach employed in Ref. 
0. In contrast to the latter, however, our present anal- 
ysis accounts for the electron dwell time in quantum 
dots and also non-perturbatively treats electron trans- 
port through the scatterers. In the limit of long dwell 
times r„±i_„ <ti 1/Tin and small channel transmissions 
/3n — > 1 (i.e. for tunnel junction arrays) our equations 
are replaced by those of Ref. 



D. Interaction correction and shot noise 

Finally, let us establish an important relation between 
interaction correction to the current and the shot noise. 
Performing summation of Eqs. H18|l with the weights i?„, 
we obtain 



1 = 



V_ 



N 

E 

n=l 



RnjCn) 

Ry. ' 



(21) 
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This formula generalizes our previous results derived for 
a coherent scattereri and a quantum dot^ to the case 
of quantum dot arrays and spatially extended disordered 
conductors. Eq. H21|) demonstrates that the interaction 
correction to the I — V curve of an array of scatterers 
scales linearly with the current noise produced by these 
scatterers. In the absence of noise the interaction cor- 
rection is identically zero, and the standard Ohm's law 
is recovered. Eq. H21|l will be extensively used in our 
subsequent calculation. 



III. RIGOROUS DERIVATION 

The phenomenological analysis presented in the previ- 
ous section clearly illustrates the relation between shot 
noise and interaction effects in electron transport. Now 
we will demonstrate that Eqs. (|9I15I19|) can also be de- 
rived within the framework of a rigorous quantum me- 
chanical procedure. This derivation will also allow to 
determine the validity range of our Langevin approach. 

We first note that for a particular case of two scatter- 
ers the above equations follow from the effective action 
analysis^ after averaging of the action over mesoscopic 
fiuctuations. Below we will see that for the case N ^ 2 
these equations yield exactly the same results as those 
derived in Ref. for chaotic quantum dots. Direct gen- 
eralization of the method^ to the case of quantum dot ar- 
rays, though technically possible, turns out to be rather 
involved since one should first establish the full quantum 
mechanical action for the whole array and then perform 
its averaging over mesoscopic fluctuations. In this case it 
appears more convenient to average the action already at 
the first stage of the calculation. In order to accomplish 
this goal wc will employ the non-linear ct— model-type of 
approach combined with the Keldysh technique. This 
method was proposed in Ref. and recently applied to 
chaotic quantum dots in Ref. 4. Below we will extend 
this technique to arrays of quantum dots. 



A. Effective action 

In the presence of electron-electron interactions 
general quantum mechanical description of both 
compact scatterersii2i2i^ii& and extended disordered 
conductorsi^iii can be formulated in terms of the effec- 
tive action which depends on the fluctuating Hubbard- 
Stratonovich fields Vi and V2 defined on the two branches 
of the Keldysh contour. In the situation considered here 
the action also depends on the fluctuating Green func- 
tion Qn which is 2 X 2 matrix in Keldysh space satisfying 
the normalization condition 



Q\ti,t2)^ dt:iQ{tiM)Q{hM)^5{ti-t2)l, (22) 



and on the fluctuating phases of the dots 



(23) 



where we defined 



- E / dt'eiVj^t') + Vj,2{t'))l2, (24) 
j=i 

=E / dt'e[V,At')~V,At')). (25) 
3=1 

Here and below ffx.y,z are the Pauli matrices in Keldysh 
space. 

The complete expression for the effective action of the 
array reads 



iS 



N-l 

f i V / dt 



e2 5„ 



2 ^ s: I ^"^n 



n=l 

N-l 



l + ^{{Qn-l,Qn}~2) 



'§i-^r.-^aA qJ(.26) 



Here the trace includes the summation over the channel 
index k while the superscript n indicates the scatterer 
number. The boundary conditions for the operators Q 



Qo{ti,t2) — — 

QN{tl,t2) 



0[tl - t2)(Tx, 



sinh7rT(ii — <2 

sinh7rr(ti — 12) 
-5{ti - t2)ax- 



(27) 



Physical observables can be evaluated by means of the 
following equation: 



{A) 



(28) 



In the FCS-type of approach^-^ the matrix Q„ has the 
form 



Qn C 



G„(CTz +icFy) - (Ta 



(29) 



where Xn is the time and space independent "counting 
field" for the n-th quantum dot. In our problem has to 
be replaced by an arbitrary fiuctuating Hermitian opera- 
tor. This observation suggests the following parametriza- 
tion of the operator Qn. 



Qn = e*^"'^^ Gn{dz +icTy)- dx e 



(30) 
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where (?„ and Wn are Hermitian operators. Here (?„ ac- 
counts for fluctuations of the electron distribution func- 
tion in the n— th quantum dot or, more generally, for 
fluctuations of the Keldysh-Usadel Green function. The 
operator Wn describes "quantum" fluctuations of the 
field Qn- It is possible to demonstrate that an arbi- 
trary operator Q„ satisfying the normalization condition 
(|22|l and being sufhciently close to the "classical" one, 

gI"'' {ffz + icTy) — cTx, can be written in the form (|30|) . We 
further note that the parametrization H30() is not identical 
to that proposed in Ref. |3 

Let us expand the action H2f)|l to the second order in 
the small operators Wn- Then we obtain 



N 



n=l 
N-1 



71=1 

N-1 



^ — Tr i^-Gn - 2 



n=l 

N 



Ot 



Wn 



E 5nTr (^-i{Gn - Gn-l)Wn + PnGn-lWnG 
(Gn-lWn) + {G 



where w„ = Wn — Wn-i- The quadratic in w„ terms can 
be decoupled with the aid of the Hubbard-Stratonovich 
transformation. One finds 



exp < ^ g„Tr ( (3nGn~lWnG 

+ [{Gn-lWnf + [GnWnf] " i^^) | 



'i:"=lTr(C„to„) 



(31) 



where we introduced the Gaussian stochastic operator C„ 
with the pair correlator 

{Cn{ti,t2)Cn{hM)) = 2g„^(t3 - t2)6{h - U) 
— dnPn [G„_i(i3, <2)G„(ti, t/C) 
+ Gn{t3, t2)Gn-l{tl,t4)^ 
- gn{l - Pn)[Gn{t3,t2)Gn{tl,U) 

+ Gn-lit3,t2)Gn-l{tl,h)]. (32) 



B. Kinetic equation and Boltzmann-Langevin 
approach 

We are now in a position to derive the equation of 
motion for the matrix Gn- In the metallic limit (?n S> 1 



it is sufficient to restrict our analysis to the least action 
condition dS/6Wn = which yields 



An 

Sn 



d_ 

dt 



i^n{t), Gn — QnGn-l + 9n+lGn+l 

- {gn + gn+l)Gn + Cn " C«+l- (33) 



Dropping the operator notations one can rewrite the 
same equation in the form 

5n \dti dt2 

+ 9n+lGn+l — {gn + 5n+l)G„ + Cn ~ Cn+l-(34) 



+ i^n{tl) - i^n{t2) Gn = gnGn-1 



Under the condition <?„ ^ 1 one can also neglect the noise 
terms ^„ in Eq. H34I) which can only contribute in higher 
orders in 1/gn- Dropping these terms one observes the 
equivalence of the equations H34|l and (|15|l . 

It is also useful to illustrate a simple relation between 
our Eq. H34II and the standard Boltzmann-Langevin 
approach^, which is frequently used, e.g., for the analysis 
of the shot noise in disordered conductors. Let us again 
define t = {ti+t2)/2 and s = ti—t2- We will assume that 
Gn varies slowly with t and is a slow function of time. 
Replacing i^n{ti) — i^n{t2) 'i^n{t)s and performing 
the Fourier transformation of Eq. (|34|l with respect to s, 
one arrives at the equation for the distribution function 

4^ fdfn dfn\ _ , ^ „ 

^ + j " •'^"^"-^ + 9n+lfn+l 

-{gn+ gn+l)fn-r]n+r]n+l, (35) 

where 2?7„ = / dse'^" C«(i + s/2, t - s/2) and 

iVnit, E)7Jrn{t', E')) = 2ngnSmnS{t - t')6{E - E') 
X {I3n[fu{t, E)hn-l{t, E) + fn-l{t, E)hn{t, E)] 
(1 - f3n)[fn{t, E)hn{t, E) + /„_i(t, E)K-l{t, S)]}.(36) 

Eqs. (|35|l . H3t)|) represent an extension of the Boltzmann- 
Langevin approach® to arrays of quantum dots. Within 
the framework of the adiabatic approximation the latter 
approach follows directly from our rigorous analysis. 



C. Excess charges 

Let us now integrate out the field $~ . Being performed 
with the action H26|l this integral gives the functional 
delta-function equivalent to the equation 5S/5'^n ~ 0- 
This equation yields 



Cgn + -7— j Vj — CnVn + Cn+lVn+1 



27re d 
Sn dt 



G„{t,t)=0. (37) 
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The last term in Eq. H37|l can be expressed via the excess 
charge g„ in the n—th dot. In the stationary case this 
charge variable is defined as follows 



1^ / dE{UE + ^,,)-ME)) 



^hnrfcj 

S„ s^o \ V 



s 



I 

ITS 



,(38) 



This definition can be applied to non-stationary situa- 
tions as well. Then Eq. P7|l is replaced by Eq. ifTTj) . 

Finally, let us consider the limit ti <2 in Eq. (|34ll . 
Setting t — {ti + 12)/2, we arrive at Eqs. (|18I19|I where 
the noise terms ^„ are identified as £,n{t) = e({t,t)/2. 
With the aid of Eq. one can check that the noise 

correlator (Cn(ii)Cn(^2)) is given by the formula @. This 
observation completes our derivation. 



IV. SHOT NOISE 

The first and immediate application of our formalism 
concerns the analysis of the shot noise in arrays of scatter- 
ers and/or quantum dots in the absence of interactions. 
We will employ Eqs. H6I12I18|I and evaluate the noise 
spectrum of quantum dot arrays in the zero frequeiic 
limit. Our procedure is similar to that applied in Ref. 
to arrays of identical chaotic cavities. 

From Eqs. (|12I13|) we obtain 



fn{E) = (1 - an)fF{E) + a^fpiE - eV), 



(39) 



where we have defined a„ — Rj/Rs- Substituting 

the result ^ into Eq. © we derive the noise spectrum 
for the n—th junction. In the zero frequency limit one 
finds 

Sn = 2(|C«|5=o) = [l ^ '^"(l ~ '^n) ~ " On-l) 

AT 

- Pn{an - an^if] — + [a„_i(l - a„_i) 
nn 

2eV eV 

+ an{l - an) + /3„(a„ - a„_i)^] — — coth — . (40) 

Un 

Finally, we note that at sufficiently low frequencies the 
term with charges (j„ in the right-hand side of Eq. H19(l 
can be neglected. In this limit the current fluctuations SI 
in the whole array are related to the current and voltage 
fluctuations across the n-th scatterer as 



(41) 



where fluctuating voltages are subject to the constraint 
En=i = 0. We obtain 51 = (l/R^) Y.^^^ i?„Cn and 



N 
n=l 



S — y R^Sn/R^. 



Then with the aid of Eq. (gUJ we get 
S 



-AT ~ 2pV pV 
(l-/3)|^+/3|^coth|l, (42) 



where (3 is the Fano factor for ID arrays of chaotic quan- 
tum dots, which is obtained in the form 



1 \ ^ i?5^ 



/?3 

1 ^Y. 



(43) 



For arrays of diffusive scatterers with (3n — 1/3 one ob- 
viously gets (3—1/3 for all values i?„. For homogeneous 
arrays with i?„ = R and /3n — f3 Eq. H43I) yields 



(44) 



This result demonstrates that in the limit — > oo an 
array of arbitrary - not necessarily diffusive - scatterers 
should behave as a diffusive conductor with f3 1/3. 
In the case of identical transmissions for all conducting 
channels Eq. 144|) reduces to that derived in Ref. |3 



V. GENERAL EXPRESSION FOR THE 
CURRENT 

Let us now turn to the calculation of the current- 
voltage characteristics in the presence of electron-electron 
interactions. According to Eq. (|^ . in order to accom- 
plish this goal it is necessary to evaluate the average value 
for the noise terms . This average would vanish identi- 
cally, {^n) = 0, should there be no dependence between 
the fluctuating voltage Vn and noise However, since 
such a dependence in general exists, the averages (^„) 
differ from zero and the interaction correction remains 
finite. 

In this paper we will restrict ourselves to the pertur- 
bation theory and provide a general expression for the 
interaction correction to the I — V curve. In what fol- 
lows we will consider noise as a small perturbation, in 
which case in the leading approximation fluctuations of 
the phase on n—th junction and the noise are related 
to each other by means of a simple formula 



S(pnit) 



1 ^ 

e 

m— 1 



drRnmit - r)^„(T). 



(45) 



An explicit expression for the function Knm{T) will be 
specifled later. Now we only point out that due to the 
causality requirement one has Knm{T < 0) = 0. Mak- 
ing use of the above relation in the lowest non- vanishing 
order one can express the average value of ^„ in the form 



iUt)) 



N 

E 

m— 1 



dti 77^d<^,„(ii) 



^ ™ / 1 <^ 



m,k—l ' 



Ut2) )(.46) 
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Here the derivative S£^n{t) / S(pm{ti) accounts for the feed- 
back of the phase fluctuations on the shot noise. FormaUy 
this effect is encoded in the Green function G„(ti,i2) 
which satisfies Eq. (|15|l and determines the noise cor- 
relator (jSJ. In the lowest non- vanishing order in ^„ it 
is sufficient to employ Eq. (|I5|I instead of H34|l . From 
the causality requirement one finds 6Cn{t)/Siprn{ti) — 
for ti > t. Utilizing this property together with a sim- 
ilar one for the function Kmkit) and making use of the 
fact that the noise variables for different scatterers are 
uncorrelated one can rewrite Eq. 146() in the form 



1 ^ f 

i^n) = J dtidt2 Kjnn{il " h) 



s{Ut)Ut2)) 



(47) 



In order to evaluate the functional derivative in Eq. 
(I47|l it is necessary to resolve Eq. (|15ll . Proceeding per- 
turbatively in Sipm and expressing the Green function 
Gn{ti,t2) in the form 

G„(ii,t2) =e-^*"(*i)+'*"(*^)f/„(ti,i2), (48) 
we arrive at the equation for the function ?7„(ti, 

+ 5„e'^"(*i)-^'^"(*^){/„_i(ti,t2) 
+.g„+ie-^'^"+i(*i)+^^"+i(*^)C/„+i(ii,t2)(49) 
With the aid of this equation we obtain 
•A ^ /■ 

m=l 

X {{Dnrn{ti - r) - A™(t2 - r))G(°L i(ti - t2) 

- {Dr^^rn-lih - t) - Z?„„„_i(t2 - r))G^^\h - ^2)1(50) 

Here G^"^ (ii — ^2) is the solution of Eq. ((TB|l obtained for 
S^n = 0. It reads 



-iT 



^.^j^^y^ (l-a« + a„e— ^^O, (51) 

where the coefficients a„ are defined after Eq. (|39|l and 
the "diffuson" £'„„j(t) satisfies the equation 



dt 



-;—i9nDn-l,m + 5ri+l-D„+l.,: 

47r ^ 



"(ffn + 9n+l)Dn^m) + (5„,„(5(t) (52) 

with the boundary conditions 

Do,m = DN,m = Dn,0 = D„^N = 0. (53) 

As before, due to causality one has Dnm{t) = for t < 0. 
With the aid of the above expressions we get 

f4^ = -.G(°)(t-t2) 

X [S{t - h) - S{t2 - ti)]e{n - m + 0) 

{ [Dn^it -h)- D„,^{t2 h)] G^^l.it - t2) 

- [Dn,m-lit - h) - D„,m~i{t2 " h)]G<^\t - t2)}(54) 



This equation in combination with Eqs. (|9I47() enables 
one to evaluate the derivative 5 {£,n{t)£,n{t2)) / 5ipm{ti) 
and derive the final expression for the current. With 
the aid of Eq. (|21|1 we obtain 



(55) 



where 51 is the interaction correction which can be split 
into two parts 51 = 5Ii + 5I2 ■ 



5Ii 



1 



N 



47r2ei?s ^ 

n,m— 1 



dtdt 



,n^T^ sineVt 



5I2 



sinh^ TTTt 

— 5n-l{am — a?i-l)-D„-l,m-l(i') 
+ 5n{am-l — an)Dn,m{t') 

- 5n{am - a„)Ai,m-i(i')] i (56) 

: —TT- y Pnlin f dtdt 7i 

7rei?| /-^ "7 sinh^TrTt 



X Kran{t " t') 
~ -C'n-l,m(i')) 



5nni5{t') " {D 

n— l,m — 1 



5n9rr 

47r 



{Dnm{t') — -Drum-l(i')) 



(57) 



Eqs. (|55I57|) represent our general result for the I ~ V 
curve of a Id array of metallic quantum dots in the pres- 
ence of interactions. One can verify that in the particular 
case of two scatterers or, equivalently, for a single chaotic 
quantum dot, Eqs. H55I57|I reduce to the expressions de- 
rived in Ref. by means of a different approach. Let 
us also note that the expression for the function Kmn (t) 
is determined by the solution of Eqs. (|17I19|) under the 
constraint H2U|) . Below we will explicitly find this solution 
for the specific case of homogeneous arrays of quantum 
dots. 



VI. HOMOGENEOUS ID ARRAY 

Consider an array formed by the scatterers and quan- 
tum dots with identical parameters. In what follows we 
set G,i = G, Cgn = Gg, gn = g, Pn ~ P, Rn — Ri 
5n = 5. In this case it is straightforward to derive the ex- 
act expressions for the functions Kmn{t) and Dnm- These 
expressions read 

2e^'^ f du e~''^* 

N I ZTT —IbJ + 

9=1 

X cos — — - cos — — — , (58) 

V N 2NJ \ N 2NJ ' ^ ' 

where we defined the impedance Zi^q 

1 



, (59) 
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and 



and 



D. 



± V / —e-'^'D sin^ 

1 



. nqn , nqm 
sm ■ 



N N 



•^9 . , l-cos(7rg/7V) 
— lUJ H — - 



(60) 



Here and below td — ^n/gS stands for the electron 
dwell time in a single quantum dot and the coefficients 
a„ reduce to a„ = n/N. Making use of the property 
Kmn{t) = Knm{t) and Dmn{t) = Dn„i{t) wc obtain 



1 + cos f dtj 

5h = -e2^ 



5=1 



2tdRN^ J 2tt 



— lm{Z^,Dl^)B{Lo,V,T) 



x-^ 1 + cos ^ /■ (ioj T 

q—l ^ 

i-(-i)%^M ^ -— ^ 

X- : tt^\/~IU}Td(2~ lUJTD), (61) 



1 + {-l)iu'^{uj) 



J 



e/5 



2N^R ^ 

q=l 



duj 
2^ 



lm{Z^qD^q)B{u,V,T), (62) 



where B{uj, V, T) = Y.±{^'^ ± ^) coth 



2T 



and u{lo) 



1 — zcj — ■\/(l — iujToY ~ 1- We observe that the second 
contribution to the interaction correction 81i scales with 
the Fano factor /3 of individual scatterers and, hence, van- 
ishes for /3 — > 0. At the same time the first contribution 
bl\ does not depend on i.e. it is universal for any type 
of scatterers. We also note that the term 6Ii differs from 
zero for all > 2 but vanishes identically in the case of 
two scatterers N = 2. 

The frequency integrals in Eqs. H61() and (|62() can be 
performed exactly with the result 



51 



2Te 

" A2 



N-l 

E 

q=l L 



/3- 



2 TUJ_ 



4 sin' 



27V J 



^ ^" ^ 2jV 
„2 _?r£ 



Sin' 111 



^sin^i^ 



4j e X - 

2sin^ 



-cos^)' 

7rp 



.7rTi?(4Csin^ II + Cg) ttTtd 2ttT ) \ ttTtd 2ttT 
1 _ (_i)P+9)sin' 

_2 TTB , _R (./-t 



2 ££ - 2 ££ 



(2 sin- 
sin' 



2lV 
2N 



R 
2td 



V ttTE (4C sin2 1^ + Cg) ttTtd 



■ (4C sin- 
27rT / \ ttTtd 



C,)(cos^-cos^)) 



2^y 



(63) 



Here we defined the function W{x) = Im[a;5'(l + x)], 
where '^(x) is the digamma function. Eq. (|63|l is the 
exact expression for the leading (in 1/g) interaction cor- 
rection to the current valid both in linear and non-linear 
in voltage regimes and for any number of scatterers TV in 
the system. 

Let us consider a physically important limit of rela- 
tively large metallic quantum dots with i?C, RCg <^ to ■ 
Making use of this inequality one can significantly sim- 
plify the general result (|53|l and find 



ryrj. N-1 

Q-1 



2sin 2lv 



1 - 



/ I sm 

V7rri?(4Csin- 

W 



■J 

2(l-(-l)") „„t2 JUL 
-1 7r(7 



1 - COS 



l' H + Cg) 

■ 2 irg 



2^, 



ieV" 
ttTtd 2ttT J 



.(64) 



In the linear in voltage regime the above expression yields 
the result for the zero bias conductance of the array G = 



1 /NR+SG. The interaction correction SG takes the form 



SG^- 



9 N-1 



^^{.ttTR (4Csin2^-f Cg 



2sin2i^ 



2(l-(--l)'') 



cot 



2 nq 
2N 



1-cos^ 



L 



2 wq 
^ 2N 



ttTt 



D 



(65) 



where we have defined L{x) — '5(1 + x) + x4''(l + x). 

Let us now briefly analyze the above results in various 
limits. The case N = 2 was already considered in details 
in Ref. @, here we will concentrate on the behavior of 
quantum dot arrays containing many scatterers ^ 1. 
In this case the expression for the interaction correction 
(|65|l can be further simplified. In the high temperature 
limit TRGrn S> 1 (where C„i =min[C, Cg]) we obtain 



SG = 



1 - 



Cn 



(66) 



6NTRC\ yCg + 4C, 
At intermediate temperatures 1 / RC, 1 / RCg ^ T ^ 



10 



1.00 



0.99 



o 

DC 
> 

DC 



0.98 



0.97 



0.96 




FIG. 2: The linear conductance G (solid lines) as a function 
of temperature T together with the differential conductance 
dl/dV (dashed lines) as a function of the applied voltage V at 
T = 0. The results are obtained from Eq. ^ for g = 1000, 
/3 = 1/3, td/RC = 10'', Cg/C = 2.5 and for three different 
numbers of barriers: N = 2,10 and 50. We identify four 
different regimes (the boundaries between them are shown 
by dotted lines): (I) saturation regime eV,T < tt'^ /2N^td, 
Eq. (1 1) diffusive regime tt'^/2N^td < eV,T < 1/td, 

Eqs. 168171)1 . (Ill) logarithmic regime of almost independent 
barriers 1/td < eV,T < 1/RC, Eqs. J6'7l'7lll and (IV) high 
temperature (classical) regime eV,T > 1/RC, Eq. Hfcifcil . 



1/td Eq. inUl yields 



5G'. 



In 



1 



2ttTRC* 



(67) 



where C* = {y^+ y/Cg+AC) /4. In the interval 
n^/2N^TD < T < 1/t£, we find 



3C(3/2) 



ttN RC* ttN A^/WTd 



e 



1 



1 
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rection to the conductance saturates, and we get 

- In l^-f 0.368 -I'ji (69) 
ttN RC* V N J TT ^ ' 

Consider now the nonlinear regime e|y| 3> T in which 
case the I — V curve is determined by Eq. (|63|) (or H64|) ). 
We will again consider the limit ^ 1 and make use 
of Eq. H64(l . At very small voltages and temperatures, 
e\V\ < Tr^/2N^TD, one finds dl/dV = 1/NR+SG, where 
SG is again defined by Eq. (|^ . At higher voltages, 
t:^/2N'^td < elP^I < l/r^i we obtain 



dV 



1 

'nr 



ttN ^ RC* 



1 



At even higher voltage, 1/td < e\V\ < l/i?C* 
differential conductance takes the form 



dl 
dV 



1 e^P 1 



(70) 



the 



(71) 



The linear conductance G as a function of tempera- 
ture and the differential conductance dl/dV at T ^ as 
a function of the applied voltage are depicted in Fig. 2 for 
different number of scatterers N in the system. One ob- 
serves that both quantities (as functions of corresponding 
variables) demonstrate a very similar behavior. In both 
cases four different regimes can be distinguished, as it 
was already specified above. Further discussion of these 
results is postponed to Sec. Vlll. 



VII. EFFECT OF ENVIRONMENT 

In order to complete our analysis let us also include the 
effect of external leads into consideration. For simplicity 
we will follow the standard procedure and assume the 
leads impedance to be purely Ohmic. This procedure 
amounts to replacing Eq. (|20|l by a more general one. 



N 



Vj + IRs = K 



(72) 



i=i 



where Vx is the voltage applied to the whole system 
"array -|-leads" , Rs is the resistance of the leads and 
/ — CVn + In is the current flowing through the leads. 
Due to current conservation the scatterer number n can 
be chosen arbitrarily here. Eq. H72|) is solved together 
with Eqs. H9I14I17I18I21|) in exactly the same way as it 
was done above in the absence of the shunt. As a result, 
the effect of Rs is accounted for by means of a simple 
replacement 



Knmit) ^ Knrn{t) + Ks{t), 



(68) where Knm{t) is defined in Eq. lf5S|l and 



where 7 ~ 0.577 is the Euler constant. Finally, in the Ks{t) 
limit of very low temperatures, T <^ 7r^/2A^^rD, the cor- 



e 

TV 



doj 
2^ 



(_,^ + 0) {-^u;+l^ + f^^ 



(73) 



(74) 
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Note that the function Ks{t) does not depend on n 
and m. After such a replacement the new expression for 
Knm{t) should be substituted into Eqs. H56I57|I and we 
arrive at the result 



I = ^ + 5I + 5Is, 



(75) 



where V = J2jLi i^j) average voltage across the ar- 

ray and 61 is defined in Eq. H63|l . The general expression 
for the additional term Sis is rather cumbersome and will 
not be presented here. Below we will only address the ef- 
fect of external leads on the linear conductance in the 
low temperature limit T < tt"^ /2N'^T£i. It turns out that 
for non-zero Rs - similarly to the case of single quantum 
dots^ - the conductance saturation is lifted and the re- 
sult (|69|l becomes incomplete. Taking into account the 
shunt contribution to the current Sis, one finds 



G 



1 



in — — 

nN RC* 



Rs 



In- 



TT Rs + NR 2N^Ttd 



where 



1 ^^^(l-(-l)«)cos2i^ 



N 
9=1 



4 Trq_ 
^ 2N 



(76) 



(77) 



The sum in (|77|l is evaluated exactly and just yields the 
Fano factor of the array H44|l . Thus, in the presence of 
an external shunt the conductance keeps decreasing log- 
arithmically with T even at very low temperatures. As 
before, this logarithmic contribution scales linearly with 
the total Fano factor of the array (3 which tends to the 
universal value 1/3 in the limit of large N. This result is 
in the agreement with our previous findingsiiS and once 
again emphasizes a direct relation between shot noise and 
interaction effects in disordered conductors. 



VIII. DISCUSSION 

In this paper we have proposed a general model which 
embraces virtually any type of disordered conductors and 
allows to account for Coulomb interaction effects in elec- 
tron transport through such conductors. Exploiting an 
intimate relation between shot noise and interaction ef- 
fects, in Sec. II and III we derived a closed set of 
Langevin-type of equations which allow to conveniently 
study electron transport in the presence of electron- 
electron interactions^*. The key idea of our approach is 
to account for modifications of the shot noise due to non- 
equilibrium effects and to self-consistently describe these 
effects and their impact on fluctuating charges and volt- 
ages inside the conductor. For the sake of definiteness 
here we focused our attention on quasi-lD conductors, 
however one can trivially extend the whole analysis to 



2D and 3D conductors as well. This generalization will 
be carried out elsewhere. 

The technique developed in this paper allows to obtain 
a general formula for the interaction correction to the 
current, Eq. H63|l . which contains all the results derived 
previously for various quasi-OD and quasi-lD disordered 
conductors and extends these results to yet unexplored 
regimes. At sufficiently high energies (exceeding the in- 
verse dwell time of a single quantum dot I/td) the scat- 
terers behave as effectively independent ones, and one 
can identify two different regimes (regimes III and IV in 
Fig. 2) described by Eqs. Ht)7l71lt)t)|) . At such energies 
the interaction correction scales with the Fano factor (3 
of individual scatterers and in a wide interval of energies 
depends logarithmically on temperature or voltage. For 
a special case of tunnel barriers — > 1 our results re- 
duce to those derived in Refs. QllJj while in the limit of 
ballistic contacts with /3 — > (or, equivalently, diffusive 
wires with point-like impurities) the interaction correc- 
tion turns out to be negligibly small in this regime. 

At energies below 1/td (regime II) scatterers located 
sufficiently close to each other become effectively cor- 
related. The number of such scatterers N^g in one 
"correlated" segment of the array grows with decreas- 
ing temperature (or voltage) as N^s ^ 1 / \JTtb (or 
A'cff ~ l/VeVTo). In this regime the system can be 
viewed as a chain of ~ N/Ncs segments, each of them 
now playing the role of a "new" independent scatterer 
with an effective conductance ^cff ~ 9/Ncs g^/TTo (or 
~ gy/eVTo)- Then the results^'M can be applied again, 
in the corresponding expression for the interaction cor- 
rection one should only substitute QcS instead of 5. In 
this case the logarithmic dependence of the interaction 
correction on temperature/voltage drops oufe^ and, e.g., 
for the linear conductance one finds SG/G ~ — /Scff/ffoff, 
where j3cS is the Fano factor of a segment with N^s 
scatterers. According to Eq. (|44|l . for sufhciently large 
TVoff S> 1 the factor /?cff approaches the universal value 
1/3, and we obtain SG/G — l/^vTro in agreement 
with the well known resul^ii and also with our rig- 
orous formula H68I) which - in addition - contains a 
temperature-independent contribution oc /3 coming from 
high energy modes. Finally, as Ncs approaches N the 
system conductance either saturates (for Rs — > 0, regime 
I) or crosses over to the low energy logarithmic regime 
(|76|l caused by additional voltage fluctuations across the 
array due to non-zero external shunt resistance Rs- 

It is also straightforward to establish a direct relation 
between the results derived here and those obtained di- 
agrammatically in the linear in voltage regimaiSiii. By 
setting N —> 00 and /? — s- 1 from Eq. H65() we reproduce 
the results^" for the interaction correction in tunnel junc- 
tion arrays, while in the limit 0—^0 the latter equation 
yields the standard resultii for diffusive wires^". The 
same equivalence can be observed at the level of gen- 
eral expressions (|61l62f) considered in the limit eV ^ T, 
iV ^ 00 and for td > RC, RGg. For /? = 1 the resultia 
follows from the sum of two terms (|61f) and Ht)2|) . while 
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for /? = the second contribution 1)62(1 vanishes identi- 
cally and the resultii is obtained only from the first term 
(|()T|l . These observations demonstrate that the reduction 
of our model to one for diffusive wires with point-like im- 
purities is achieved by setting /? ^ 0. In the latter case 
N coincides with the total number of impurities in the 
wire. 

At last, let us briefly summarize the applicability con- 
ditions for our results. As it was already discussed above, 
our Langevin approach is justified in the metallic limit 
(7„ 3> 1. Under this condition our technique should ac- 
count for all essential processes except for subtle instan- 
ton effects which may show up only at exponentially 
low energies. An obvious necessary (though possibly 
not sufficient) validity condition of our results derived 
in the linear in voltage regime is 6G/G <C 1. While at 
high enough temperatures this inequality is automati- 
cally fulfilled in the metallic limit gn ^ 1, at the lowest 



energies/temperatures a much more stringent condition 
5s — 27r/e^i?i; » 1 has to be satisfied. The latter condi- 
tion is inevitably violated for large number of scatterers 
N in which case a non-perturbative analysis becomes nec- 
essary in the low energy limit. This analysis is beyond 
the frames of the present paper. In the non-linear regime 
and at sufficiently high voltages the applicability range 
of our results can be additionally restricted by electron 
heating effects which we do not address in this work. 
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In order to avoid misunderstandings we would like to em- 
phasize again that our present approach accounts for the 
interaction correction to the I—V curve and is not meant to 
describe weak localization effects and interaction-induced 
dephasing^''^'^^ which we do not address in this work. 
In this case the combination Tt]^ occurs under the loga- 
rithm, where is the total dwell time for Nbb scatterers. 
This combination reduces to a temperature independent 
constant of order unity since r|f ~ l/T. 
In order to observe the equivalence between our Eq. 16511 
in the limit (3 —* and the results— for the interaction cor- 
rection in quasi- ID diffusive wires it suffices to introduce 
the diffusion coefficient D = a? /2td, where a is the size of 
a single quantum dot. 



